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The dynamic atomic polarizability describes the response of the atom to incoming electromag¬ 
netic radiation. The functional form of the imaginary part of the polarizability for small driving 
frequencies uj has been a matter of long-standing discussion, with both a linear dependence and 
an Lo^ dependence being presented as candidate formulas. The imaginary part of the polarizability 
enters the expressions of a number of fundamental physical processes which involve the thermal 
dissipation of energy, such as blackbody friction, and non-contact friction. Here, we solve the long¬ 
standing problem by calculating the imaginary part of the polarizability in both the length {“d-E”) 
as well as the velocity-gauge (“p ■ A ") form of the dipole interaction, verify the gauge invariance, 
and hnd general expressions applicable to atomic theory; the ui^ form is obtained in both gauges. 
The seagull term in the velocity gauge is found to be crucial in establishing gauge invariance. 

PACS numbers: 31.30.jh, 12.20.Ds, 31.30.J-, 31.15.-p 


I. INTRODUCTION 

The “susceptibility” of an atom toward the generation 
of an atomic dipole moment is described by the dynamic 
atomic (dipole) polarizability. Alternatively, the atomic 
polarizability describes photon absorption from a field of 
photons and subsequent emission of a photon into the 
same or a vacuum mode of different wave vector but the 
same frequency. The external oscillating field can be a 
laser field (and this is the intuitive picture we shall use 
in this paper). In this case, the real part of the polariz¬ 
ability describes the ac Stark shift [l|. The atomic po¬ 
larizability has both a real as well as an imaginary part. 
Physically, this can be illustrated by considering the rel¬ 
ative permittivity er(w) of a dilute gas and its relation to 
the dynamic dipole polarizability a(a;) of the gas atoms, 

er{u}) = 1 +—a{uj), (1) 

eo 

where eo is the vacuum permittivity and Ny is the vol¬ 
ume density of gas atoms 0. The Kramers-Kronig re¬ 
lations dictate that the real part of e(a;) should be an 
even function of the driving frequency w, while the imag¬ 
inary part of e(a;) cannot vanish and should be an odd 
function of w. These considerations are valid upon an 
interpretation of the dielectric constant in terms of the 
retarded Green function Gr which describes the relation 
of the dielectric displacement D{r,t) to the electric field 
E{r,t), 

poo 

D(f,t) = eoE(f,t)eo / dr Gr{t) E{r,t - r). (2) 
Jo 

The Fourier transform is 


literature The conventions used here and also the 

result of our calculation of the non-resonant contribution 
to the imaginary part, to be reported in the following, are 
consistent with the interpretation as a retarded Green 
function and, notably, with the conclusions of Ref. Q 
[see Eq. (2) and the text following Eq. (31) of Ref. 0]- 
The imaginary part of the polarizability enters the 
description of processes such as the black-body fric¬ 
tion [ 3 , @ , where an atom is decelerated by interaction 
with a thermal bath of black-body photons. The black- 
body (BB) friction force Ebb = —"Obb v is linear in the 
velocity. 


00 
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Here, [3 = l/{kBT) is the Boltzmann factor, c is the 
speed of light, and h is the quantum unit of action (re¬ 
duced Planck constant), while the integration is carried 
out over the entire blackbody spectrum of angular fre¬ 
quencies Lo. This force is due to dissipative processes; 
the atom absorbs an incoming blue-shifted photon from 
the front, and emits photons in all directions. The fric¬ 
tion force at a distance Z above a surface composed of 
a dielectric material, due to the dragging of the mirror 
charges inside the electric with frequency-dependent di¬ 
electric constant e(a;), is calculated as Eqf = —rjQF v. 
According to Ref. [1^, the result for the non-contact 
quantum friction (QF) coefficient in the relation is given 
by a Green-Kubo formula, 
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Gij(w) =- a{uj ). (3) 

eo 

The symmetry properties of the imaginary part of the 
polarizability have been discussed at various places in the 


The formulas (|3]) and (O raise the question what could 
be deemed to be the most consistent physical picture be¬ 
hind the imaginary part of the polarizability. The imag¬ 
inary part of the polarizability involves a spontaneous 
photon emission process, and this spontaneous emission 
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can only be understood if one quantizes the electromag¬ 
netic field. From a quantum theory point of view, the 
reference state in the calculation of the leading-order po¬ 
larizability is given by the atom in the reference ground 
state and ul photons in the laser field. The virtual states 
for the calculation of the polarizability involve the atom 
in a virtual excited state, and n/, ± 1 photons in the laser 
field. The dipole coupling leads to the creation or anni¬ 
hilation of a virtual laser photon. 

Let us try to calculate the imaginary part of the po¬ 
larizability on the basis of an energy shift calculation. 
In many cases, the imaginary part of an energy shift de¬ 
scribes a decay process, with the spontaneous emission of 
radiation quanta. If one includes the energy of the emit¬ 
ted quanta into the energy balance, one sees that the 
final quantum state of the decay process actually has the 
same energy as the initial state. We consider as an ex¬ 
ample the spontaneous emission of a photon by an atom. 
The initial quantum state has the atom in an excited 
state, and zero photons in the radiation field. In the fi¬ 
nal state of the process, the atom is in the ground state, 
and we have one photon in the radiation field, with an 
energy that exactly compensates the energy loss of the 
atom. The imaginary part of the self-energy of an ex¬ 
cited state describes the decay width of the excited ref¬ 
erence state, against one-photon decay m- This consid¬ 
eration has recently been generalized to the two-photon 
self-energy [I^, which helped clarify the role of cascade 
contributions which need to be separated from the coher¬ 
ent two-photon correction to the decay rate [l^ . 

How can this intuitive picture be applied to the imag¬ 
inary part of the polarizability? From a quantum theory 
point of view [l|, the initial state/reference state in the 
calculation of the polarizability involves the atom in the 
ground state, and ul photons in the radiation field. The 
laser photons are not in resonance with any atomic tran¬ 
sition, but rather, in a typical calculation, of very low 
frequency. One may ask what is the energetically degen¬ 
erate state to which the reference state could “decay”. 
In the case of the ground-state atomic polarizability, one 
may additionally point out that the atom already is in 
the ground state and cannot go energetically lower. How 
could the “decay rate” be formulated under these con¬ 
ditions? The answer is that a quantum state with the 
atom in the ground state, 71^ — 1 laser photons of energy 
loli and one photon of wave vector k and polarization 
A, with energy = ujl (but not the same polariza¬ 
tion or propagation direction) is energetically degenerate 
with respect to the reference state (with the atom still in 
the ground state, and ul laser photons). The reference 
state in our calculation will thus be the product state of 
atom in the reference state (in general the ground state) 
\4>), the laser field in the state with ul photons, and zero 
photons in other modes of the quantized electromagnetic 
field. We denote this state as 

|0o) = . (6) 

The calculation of the energy shift then proceeds in anal¬ 


ogy to the self-energy calculation: One inserts radiative 
loops into the diagrams that describe the ac Stark shift 
and calculates the decay to the state 

\(j)f) = \(t),nL-l,lkx) ■ (V 

The notation adopted here involves the occupation num¬ 
bers of the laser mode (subscript L), and the mode with 
wave vector k and polarization A. 

The remainder of this paper is organized as follows. 
In Sec. El we discuss the derivation of the imaginary 
part of the polarizability in the length gauge, where the 
interaction of the atom with the laser field, and also 
the interaction of the atom with the quantized radia¬ 
tion field, are modeled on the basis of the dipole inter¬ 
action —e r • E. Gauge invariance (see Sec. ED is used 
as a method to verify our results. In the velocity gauge, 
we use the “dipole” coupling —ep-A/m and the “seagull 
term” e‘^A‘^l{2m). Here, A is the vector potential, while 
E = —dtA is the electric field. We employ the dipole ap¬ 
proximation throughout the paper and work in SI mksA 
units (with the exception of a few illustrative unit con¬ 
versions in Sec. IlYl). 


II. LENGTH GAUGE 


We start with the length-gauge calculation. The laser 
is z polarized, and we consider the electric dipole interac¬ 
tion which is the dominant interaction for an atom. The 
electric field operators for the laser field (subscript L) 
and for the quantized field interaction (subscript I) read 
as follows. 



Hl= —czEl^ Hj = —ef-Ei. (8c) 


We employ the dipole approximation, setting the spatial 
phase factors exp(ifc • r) equal to unity. The use of an ex¬ 
plicit normalization volume Vl for the laser mode (and 
V for the modes of the quantized field) allows for a con¬ 
sistent normalization of the laser intensity [see Eq. m 
below]. The sum over the modes of the vacuum in Eli 
explicitly excludes the laser mode, in the sense of the 
requirement k X ^ L. However, the absence of the one 
mode does not enter the matching condition 


E = v 


d^fc 

(^’ 


(9) 


because the one highly occupied laser mode does not fea¬ 
ture prominently in the sum over all modes, and the sum 
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over polarizations can be carried out using the relation 


k\ k\ 

\ 




( 10 ) 


The quantity Vl is the normalization volume for the laser 
field, whereas V is the corresponding term for the quan¬ 
tized field mode. The laser field intensity is 


absorbed into a renormalized reference state energy E.) 
In the matrix elements used in this paper, the reference 
state is the full quantized-field reference state \(j)o) un¬ 
less indicated differently. After tracing out the photon 
degrees of freedom, one obtains 
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The unperturbed Hamiltonian is given as follows. 


Ho = 

Ha + Hem j 

(12a) 

Ha = 

^ ^ Em l^m) (0m | : 

m 

(12b) 

Hem = 

I] a4a “fe a + a+ gl . 

k\^L 

(12c) 


Here, the subscript A stands for atom, while the subscript 
EM denotes the (quantized) electromagnetic field. The 
unperturbed state is given by 

|0o) = If/),'reL,0), Ha\ 4>) = E\(l)). (13) 

denoting the atom in the reference state \(j)) (notably, 
the ground state), ul photons in the laser mode, and the 
vacuum state for the remaining modes of the electromag¬ 
netic field. We should add that the Fock state of the laser 
field is just a calculational device for our derivation of the 
imaginary part of the polarizability. (Of course, in a Fock 
state, the occupation number is precisely ul, which is in 
contrast a coherent state which is a superposition of Fock 
states.) 

The atomic Hamiltonian (Schrddinger Hamiltonian) is 
denoted as Ha^ and we do all calculations for hydrogen, 
here, whereas for many-electron atoms, in the nonrel- 
ativistic approximation, one would replace ^ —>■ Za, 
where a denotes the oth electron. All results below will 
be reported for atomic hydrogen with a single electron 
coordinate, but the results hold more generally, with the 
only necessary modification for a many-electron atom be¬ 
ing the addition of the coordinates of the other electrons. 
The reference state fulfills the relations 


which is a matrix element that involves only the atomic 
Green function, and is evaluated on the atomic reference 
state \(j)). With the identification m of the laser inten¬ 
sity, the second-order ac Stark shift can be expressed in 
terms of the dipole polarizability a(w), which we write 
as 




a{ui) 
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(17a) 


1 


E i Swl 



(17b) 


for a radially symmetric atomic reference state \(j)) which, 
in general, is the ground state. The a;* are the Cartesian 
components of the electron coordinate (the summation 
convention is used for i), while the sum over the two 
terms with both signs Eoj will be encountered several 
times in the following; it is denoted by the symbol 
here. 

The polarizability of an atom can be written as the sum 
of “positive-frequency” and “negative-frequency” compo¬ 
nents a+{uji,) and a-{uji,), according to the formulas 


a{iOL) 

a+{ujL) 


a-{ujL) 


a+{uJL) + a-iuJL), 


(18a) 
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Ha — E + hu)L + ie 

1 

Ha — E — hjjjL — \ e 


(18b) 


(18c) 


Ho\(j)o) = Eo\(j)o), Eo = E + nLhujL, (14) 

where is the number of laser photons. 

In the quantized formalism [l| , the second-order energy 
shift, which we use for normalization, can be calculated 
as 

where the prime denotes the reduced Green function. 
(Replacement of 77^ by Hj leads to the low-energy part 
of the Lamb shift, see Ref. M- Here, the Lamb shift is 


where the infinitesimal imaginary parts are introduced in 
accordance with the paradigm that the poles of the po¬ 
larizability, as a function of ojl, have to be located in the 
lower half of the complex plane. This is consistent with 
the fact that the polarizability corresponds to a causal 
(retarded) Green function. We take note of the identity 

—— = (P.V.)—-f i7r^(a:). (18d) 

X — le X 

We the help of this relation, one easily shows that the 
resonant, tree-level imaginary part of the polarizabilty 
Im[aij(a;L)] is an odd function of the driving frequency 
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Ini[aij(a;L)] = ^ |a;*| 5(£'„ - E - huL) 


- X! K'^ 1^1 - E + fu^L) 

n,i 

TT fnO S{En — E — HlOl) , . 

= o -;?-^ ^ 


— £/ 


(18e) 



a 



b 
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where E^ — -E is the resonant frequency for the excita¬ 
tion of the ground state to the excited virtual state |(/)„) 
of energy E„, and we have implicitly defined the dipole 
oscillator strength /„o according to Ref. [T^ . The second- 
order contribution to the polarizability Im [a^iujL)] is in 
itself gauge invariant, and it constitutes a sum of resonant 
peaks, which are relevant as the laser frequency is tuned 
across an atomic resonance. Note that the contribution 
Im[Q;^(a;L)] corresponds to the cutting of the following 
Feynman diagram, 



according to the Cutkosky rules M- Wiggly lines denote 
photons (here, the interactions with the photons of the 
laser field). 

However, there is an additional contribution to the 
imaginary part of the polarizability, created by the inser¬ 
tion of a virtual photon into the second-order diagram, 
as shown in Fig. [T] The terms with paired interactions 
with the laser and the radiation field, in the fourth order, 
are given by 

= - {Hi G'(Eo) Hl G'(Eo) Hl G'(Eo) Hi) 

- {Hl G'(Eo) Hi G'(Eo) Hi G'(Eo) Hl) 

- 2 {Hi G'(Eo) Hl G'{Eo) Hi G'{Eo)Hl) , (19) 

where G'{Eq) = {[/{Hq — Eq)'] is the reduced Green func¬ 
tion. Written in this form, the terms correspond to the 
entries in Fig.[Tl (a) and (b), and the last term to the sum 
of (c) and (d), while the photon loop involves a photon in 
the mode k X. The propagator denominators in Fig. [TJa) 
in the “outer” legs of the electron line read as H—E+loj:^ 
(twice) because the spontaneously emitted photon is part 
of the virtual state of atom-ffield. By contrast, the prop¬ 
agator denominators in the “outer” lines in Fig.[T](b) read 
as H — E — ujl (twice) because the spontaneously emitted 
photon is not present. In all cases, one photon is “taken 
from” the laser field. For Figs. [IJc) and (d), we have a 
mixed configuration with two propagators, with denom¬ 
inators {H — E + uj) and {H — E — oj). 



c 

FIG. 1: (Color online.) Feynman diagrams contributing to the 
imaginary part of the polarizability in length gauge. Wiggly 
lines denote photons, straight lines denote the bound electron. 
A laser photon (interactions with the external laser field are 
denoted with a cross) is absorbed from the laser field, while 
the self-energy loop describes the self-interaction of the atomic 
electrons. When the virtual state denoted by the internal line 
becomes resonant with the reference state, a pole is generated 
in the integration over the degrees of freedom of the virtual 
photon (Cutkosky rules [T^ , indicated by the vertical dashed 
lines). 


In the spirit of Ref. [13, the diagram (a) in Fig. [1] 
entails the following energy shift. 
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00 ) . (20) 


We now use the matching condition (0 and calculate the 
sum over polarizations using Eq. m- Isolating the terms 
which describe the absorption of a photon from the laser 
field and emission into the field mode k A, we obtain the 
following expression, 


SEa-^ -e^ ^ 
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d^k nifujJL 
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( 21 ) 
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The ^ sign implies that the terms which involve the vir¬ 
tual state \(j),nL — still need to be isolated. We 

have traced out the photon degrees of freedom. We can 
isolate the imaginary contribution created by the virtual 
state 10, n-L — 1,1^:;,^), with = ujl, as follows, 

1 

Ha - E - hioL + 

1 

Ha - E - hujL + - ie 


reads as 
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( 22 ) 


Here, an infinitesimal imaginary part has been intro¬ 
duced according to the Feynman prescription for the 
positive-energy, virtual atomic states. Keeping track of 
the prefactors, one finally obtains 
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Ha — E + fujjL 



where the Einstein summation convention is used for re¬ 
peated indices. The term in curly brackets, upon com¬ 
parison with Eq. (II3, can be identified as the imaginary 
part of the polarizability associated with diagram (a) of 
Fig. dJ The imaginary part Im((5Ea) of the energy shift 
associated with Fig. [TKa) is negative, as it should be for 
a decay process, while the imaginary part of the polariz¬ 
ability is positive [see the overall minus sign in Eq. (El)]. 

Diagram (b) in Fig. [1] generates an imaginary part in 
much the same way, but the “outer” virtual states have 
one virtual laser photon less than the intermediate state 
which generates the imaginary part, hence 


iliji{dEi,) 
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Finally, the sum of diagrams (c) and (d) generates the 
following expression. 


After matching with Eq. El, it can be summarized in 
the following, compact result, 

Im[a(wL)] = lm.[aR{uJL)] + -—^ [a(u)L)Y > (27a) 

DTreoC^ 


where Im[Q:/{(a;i)] is the sum of the resonance peaks, ac¬ 
cording to Eq. (IlSel) . and in the last term, the square of 
the polarizability is calculated using propagators without 
any added “damping terms” for the virtual state ener¬ 
gies, i.e., under the identification [a{oJL)]‘^ = Re[Q;(wL)]^ 
[see Eq. (l26l) ]. The result (l27l) is manifestly odd in the 
argument lol-, as it should be [see Eq. ©]• 

Let us dwell on the precise form of this result a little 
more. As shown below in Eq. (I48L the second term on 
the right-hand side of Eq. (I27a is a correction of rela¬ 
tive order where q;qed is the fine-structure con¬ 

stant. Eurthermore, the contribution given in Eq. (IlSel) 
can be identified, in terms of a quantum electrodynamic 
formalism, as the “tree-level” contribution to the imag¬ 
inary part, or, as the imaginary part of the “tree-level” 
polarizability itself. If we absorb in our defini¬ 

tion of aTi^{uj) relativistic corrections to the dipole po¬ 
larizability of relative order q;qeq (see Ref. [3) and all 
radiative corrections without cuts in the self-ene rgy ph o- 
ton lines (of relative order aggp, see Refs, [l^ [l9lEI| l. 
then we can write Eq. (I27al) alternatively as 


Im[a(a;)] Im[Q;TL(i^)] 


a{ui) 


aTL(w) 


67reoc^ 


q;tl(w) + Cl(aQED) • 
(27b) 


We can “sum” the second term into a denominator. 


Im[a(a;)] « Im 


1 




(27c) 
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(25) 


The end result for the imaginary part of the polarizabil¬ 
ity, generated by the fourth-order diagrams in Fig. [TJ 


The latter functional form is in agreement with various 
results for polarizabilities of nanoparticles (not atoms) 
and other structures, corrected for “radiation damping” 
[see ^ (5) of Ref. ^ (27) of Ref. [H, Eq. (^of 
Ref. [H, Eq. (4) of Ref. [H, and Eq. (34) of Ref. [3]. 
The term from the denominator in Eq. (l27cD has the 
required functional from for “radiation damping” due 
to the Abraham-Lorentz radiation reaction force [see 
Eq. (1) and (10) of Ref. [l^, and Refs. [11,123]. The tree- 
level term is identified with the “damping by absorption” 
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of the nanoparticle, 


Im[aTL(w)] « Im[a(a;)] - [aTL(w)]^ . (28) 

or, as the “non-radiative” imaginary part in the sense of 
Ref. [ 2 ^. Note that the “tree-level” for a nanoparticle 
term is denoted as ao(<^) in Ref. [l^, which however, for 
an atom, would be associated with the “monopole” polar¬ 
izability which is vanishing (canonically, one denotes by 
ai{uj) the 2^-pole polarizability of an atom, see Ref. [ 13 ). 

A possible interconnection of the imaginary part of the 
polarizability with the square of t he p olarizability has 
previously been formulated in Refs. [13 and [3l|. A term 
with the square of the polarizability is added to the imag¬ 
inary part of the polarizability in Eqs . (G2) and (G3) of 
Ref. [23 [see also Eq. (49) of Ref. [13. Within QED, the 
term with the square of the polarizability is identified as 
being due to a one-loop correction, while the resonant 
term is added to the results given in Refs. [H, [H, |3l| 
and constitutes a tree-level contribution. 


III. VELOCITY GAUGE 


The question of the choice of gauge in electrody¬ 
namic interactions has been raised over a number of 
decades [13, with a particularly interesting analysis be¬ 
ing presented in Ref. [13 • In general, and in accor¬ 
dance with the now famous remark by Lamb on p. 268 
of Ref. 3J|, the wave function of a bound state there¬ 
fore preserves its probabilistic interpretation only in the 
length gauge, and it is this gauge which should be used 
for off-resonance excitation processes [13113 • According 
to Table 1 of of Ref. [13 , the dipole matrix elements for 
dipole transitions from state |i) to state [/) in the dipole 
approximation are related by the 
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—p ■ A 

m 


huj 




er • E 


(29) 


as can be shown with the help of the commutation rela¬ 
tion P = ^ [Ha , r]. The two forms are equivalent only at 
resonance hu = Ef — Ei. In quantum dynamic calcula¬ 
tions, which often involve nonresonant driving of atomic 
transitions, according to Lamb’s remark on p. 268 of 
Ref. [13 , “a closer examination shows that the usual 
interpretation of probability amplitudes is valid only in 
the former [the length] gauge, and no additional factor 
[{Ef — Ei)/{hjjj)] occurs.” Namely, the momentum oper¬ 
ator p = —\hV describes the mechanical momentum of 
the matter wave in the length gauge, where the kinetic 
and canonical momentum operators assume the same 
form [13113 • In thfi velocity gauge, the canonical momen¬ 
tum p — e A assumes the role of the conjugate momentum 
of the position operator in the equations of motion. One 
can also argue that the electric field used in the length- 
gauge interaction is gauge invariant, while the vector po¬ 
tential in the velocity-gauge term is not [illll,[13S3- 


However, for processes such as energy shifts, where the 
initial and final states have the same energy and the sum 
of the exchanged photon energies need to add up to zero, 
both gauges should lead to the same results (see Sec. 3.3 
of Ref. [ 13 , and Refs. [Hl-Q). While limits of gauge 
invariance have recently been discussed in Ref. [43, it 
is a nontrivial exercise to check the gauge invariance of 
the imaginary part of the polarizability. Similar questions 
have been discussed in the context of the QED radiative 
correction to laser-dressed states (Mollow spectrum, see 
Refs. [Il,[li-|43). 

In the velocity gauge, the interactions are formulated 
in terms of the vector potentials A^ and A/, which de¬ 
note the laser field and the interaction with the quantized 
radiation field. We denote the corresponding interaction 
Hamiltonians, in the dipole approximation, by the low¬ 
ercase letters and hi and write 



(30a) 

(30b) 

(30c) 


From the seagull term, we isolate the terms which involve 
at least one interaction with the laser mode and write 


u 2 {-^L + AjY ^ ^ 

hs = e -r- > hLi + h^L , 


hii = e 


2m 

2 Al ■ Ai 

m 




In the length gauge, the second-order shift in the laser 
field is given by Eq. m- In the velocity gauge, one ob¬ 
tains an additional term, as follows. 


= (hLi) 



1 

Ho-Eo 



(32) 


We trace out the laser photon degrees of freedom and 
obtain 
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Ha - E - hwL 


El 

m 



(33) 


where Pz is the z component of the momentum operator. 
One might conclude that this form is manifestly different 
from Eq. m- Using the operator identity 
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i 
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[H - E + hwL,r''[ , 


(34) 
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a 
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c 
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From the diagram in Fig.[2^a), we have the energy shift 


6£a = - ( 4>o 


hLi 


Hn-Er 


■ hLi 


(36) 


After tracing out the photon degrees of freedom, the 
imaginary part can be extracted as follows, 


lin{S£a) = -i 


II 

2eoc 67reo irfitOL 


(37) 


From the diagram in Fig.[^b) and (c), we have the energy 
shift 


5£b+c = 2 {(j)Q 


hi 


Hr. -El 


■ hL 


Hn-Er 


■ hLi 


</>0 


(38) 


The imaginary part of the energy, generated by the dia¬ 
grams Fig. [2](b) and (c), can conveniently be written as 
follows. 


Im((5fb+c) = 


II 

2eoC 37reoc^m 


WL 


(39) 



1 

p* 

m H — 

E + hwL 
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Finally, from the diagrams in Fig. EJd) and (e), we have 
the energy shift 


FIG. 2: (Color online.) Additional Feynman diagrams involv¬ 
ing the seagull graph which have to be considered in the ve¬ 
locity gauge, in order to show the gauge invariance of the 
imaginary part of the polarizability. 


one may show the relation 


1 

3 





p’- 1 p* 

m H — E + hujL rn 



H — E + hujL 



where r ^ = a;* x*. With the help of this relation, and 
assuming a spherically symmetric atomic reference state 
|(/)), one may confirm that the second-order shifts (1551) 
and (13 are equal. 

In the additional diagrams which persist in the velocity 
gauge (see Fig. El), the seagull graph enters “in disguise”, 
with the virtual photon of the self-energy loop and the 
laser photon emerging from the same vertex. We may 
anticipate that diagram (a) in Fig.[2]generates a constant 
term, while diagrams (b) and (c) generate terms with 
one propagator, with a propagator denominator of the 
form H — E + ojl (the spontaneously emitted photon is 
present in the “outer” lines). By contrast, the diagrams 
in Fig. [5] (d) and (e) generate propagator denominators 
H — E — ojl for the internal line of the diagram. 


5£d+e — 2 ( (po 


hi- 


Ho — Er\ Hr, — El 


■ hLi 


Po) (40) 


where the sequence of the dipole coupling and hj is re¬ 
versed in comparison to the diagrams in Fig.[2](b) and (c). 
The imaginary part of the energy shift is the same up to 
a sign change in the frequency term in the propagator 
denominator. 


lu\{5£d+e) = i 
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(41) 



The diagrams in Fig.[2]all involve at least one occurrence 
of the combined laser-quantized-field seagull term 

Finally, we can write the additional imaginary part 
due to the seagull terms, denoted by Im(df 2 ) = 
Im(dfo+f,+c+d+e), in the velocity gauge as follows. 


lm{5£2) = — i 
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We have used the identity (1351) . To this result we have to 
add the contribution of the diagrams in Fig. [U this time 
evaluated with the dipole Hamiltonians replaced by their 
velocity-gauge counterparts and Hj hj), 


ilm(5fi) 


. II 1 1 

2eoc hTTCoc^ ujl 



p* 1 p* 


3 

TO El — E ± kuJL TO 

vj 


(43) 


Again, with the help of (I55|) . we can write this as 

ilm(^fi) = - ^ 3 - \-2ujLa{uj)\ 

Zcoc OTTeoC'^m iitilul J 

Finally, one obtains in the velocity gauge 

Im(Jfi) + Im(5£:2) = a;i [a{u:i.)f , (45) 

confirming the result ((26ll once more [the oscillator 
strength in Eq. (TTSl) is in itself gauge invariant, and the 
resonant contribution (IlSep needs to be added to the off- 
resonant contribution (HSl) ]. 

IV. CONCLUSIONS 

Our main result concerns the imaginary part of the 
atomic polarizability, which is the sum of a resonant term 
Im[a/{(a;L)] and of an off-resonant driving of an atom 
with the dynamic dipole polarizability a{uj). According 
to Eq. (I27l) . the result is given as follows, 

Im[a(wL)] = lm.[aH{ujL)] + 7 —^ Re[a(wL)]^ . (46) 

bTreoC'^ 

The clarification of the functional dependence of the non¬ 
resonant contribution to the imaginary part for small 
to has been a matter of discussion in the past (see 
Chap. XXI of Ref. [i^ and Refs. [^, E^), with both 
a linear dependence on ui and an ui^ dependence being 
discussed as candidate formulas. Initially, one would as¬ 
sume that the uj^ dependence is favored in the “length 
gauge”, whereas the linear term is obtained in the “ve¬ 
locity gauge” (see the discussion in Ref. 0), but the 
formulation presented here in terms of the imaginary 
part of a (necessarily gauge-invariant) fourth-order en¬ 
ergy shift removes the ambiguity. So, our result (H51) set¬ 
tles the question. While the calculations have been de¬ 
scribed for a single-electron atom (hydrogen atom), the 
results hold more generally because one may simply add, 
in the calculation of the polarizability, the dipole cou¬ 
plings of the other electrons, according to the replace¬ 
ment —>■ xl^ where a denotes the subscript of the 


electron. The full inclusion of all diagrams given in Figs.[T] 
and [5] is crucial in confirming the gauge invariance (see 

Secs. El and uni. 

Our results indicate that the imaginary part of the 
atomic polarizability cannot be obtained on the basis of 
a replacement in the propagator denominators of the po¬ 
larizability function, 


a{ijj) = 




3 Em — EE tujj 

|{0|^ |0m)| 

3 ^ Em — ^Tm — E ±fux 


(47) 


Beyond tree level, the latter prescription cannot possi¬ 
bly lead to a consistent result for the imaginary part of 
the polarizability, no matter how one extrapolates the 
decay width F^ = Fm(aj) off resonance. However, we 
may at least observe that if the Fm represent the one- 
photon decay widths of the virtual states, then the re¬ 
placement gll) generates the right order-of-magnitude 
for the imaginary part off resonance, in agreement with 
Eq. (H51) . 

The specification of the real part of the square of the 
polarizability in Eq. (1461) serves to maintain the rele¬ 
vant symmetry: Namely, Im[a(wL)] needs to be an odd 
function under a sign change of its argument. However, 
the off-resonant (one-loop) contribution to the imaginary 
part of the polarizability is in itself a radiative correc¬ 
tion and should thus be suppressed by powers of the 
fine-structure constant aqED- Hence, to leading order 
in aQED, we can leave out the additional specification 
of the real part in the second term in Eq. (H51) . To put 
this statement into context, we switch to natural units 
with = c = eo = 1 and refer to Ref. [i^ where it 
is shown that the imaginary part of the one-loop self¬ 
energy describes an imaginary energy contribution of or¬ 
der oiQgQ TO, where aqED is the fine-structure constant 
[in contrast to a{uj) which denotes the polarizability]. 
The radiative correction to the energy is of relative order 
ofQED in comparison to the Schrodinger energy; the latter 
is of order agg]-) in natural units. 

Let us apply this program to the polarizability. Taking 
into account that the Bohr radius is of order I / (aqED xn) 
in natural units (i.e., about 137 times larger than the 
reduced Compton wavelength of the electron), we obtain 
the following order-of-magnitude estimates: 

I X aQED 1 1 

a(w) ~ 7 -^- = 7 -DT . (48) 

(aqEDw)^ “qED”^ (aqED^i)'^ 

In order to perform the estimate for Im[a(a;)] ^ to“^, we 
need to take into account that a typical atomic driving 
frequency is of order ujl ~ ctqED ^nd the ratio again 
is of the order of 


Im[a(a;)] 

a{ui) 


3 

O^qED I 


( 49 ) 
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as it should be for a one-loop radiative correction [here, 
Im[a(w)] is restricted to the second term in Eq. (H^ ]. 

Because of the general nature of the obtained result, 
it is useful to discuss the conversion to other unit sys¬ 
tems. In natural units (n.u.), with h = c = eo = the 
result (|46]) simplifies to 


Im[a(w)]|„ u. = Im[aji(a;)]|„ -b ^ . 

(50) 

Using the relation = 47raQEDj we can write the dipole 
polarizability, in natural units, as follows. 


q;(w)| 


47rQ!QED 

“qed 


X 



aQED m x-' — -—-aQED m x- 

H — h/ zh UJi, 



Here, the sum over the signs ± denotes the two terms 
generated by virtual absorption and emission of the laser 
photon, as in Eq. (I17b|) . For an atom, all quantities 
in the term in curly brackets are of order unity; e.g., 
OiQETi'nix^ = x^jaQ = fP where oq is the Bohr ra¬ 
dius. The numerical values (reduced quantities) associ¬ 
ated with the polarizability, in atomic units (a.u.) and 
natural units (n.u.), are thus related by 

«(‘^)ln.u. = «(‘^)la.u. • (52) 

^QED 


Atomic transition frequencies are measured in Hartrees 
in atomic units, and thus we have 

‘^iln.u. =aQEDW ■ (53) 

So, in atomic units, our main result given in Eq. (EZD 
and (I50|) reads as 


2q,3 

Im[a(w)]la.u. = Im[aij(‘^)]la.u. +-“^(‘^)]|au ■ 

(«) 

The imaginary part of the polarizability enters the de¬ 
scription of a number of dissipative processes, such as 
the quantum friction due to interaction with a bath of 
black-body photons dl]), or, with a dielectric surface [see 
Eq. ([5])] . The imaginary part of the polarizability is man¬ 
ifestly nonvanishing off resonance and, for small driving 
frequency w, is proportional to 


ACKNOWLEDGMENTS 

Helpful conversations with M. De Kieviet are grate¬ 
fully acknowledged. The authors (U.D.J.) wish to ac¬ 
knowledge support from the National Science Foun¬ 
dation (Grants PHY-1068547 and PHY-1403973) and 
(K.P.) from the Polish National Science Center (NCN, 
Grant 2012/04/A/ST2/00105). Early stages of this re¬ 
search have also been supported by the Deutsche 
Forschungsmeinschaft (DEG, contract Je285/3-2). 


[1] M. Haas, U. D. Jentschura, and C. H. Keitel, Am. J. 
Phys. 74 , 77 (2006). 

[2] J. D. Jackson, Classical Electrodynamics, 3 ed. (J. Wiley 
& Sons, New York, NY, 1998). 

[3] D. L. Andrews, L. C. Davila Romero, and G. E. Stedman, 
Phys. Rev. A 67 , 055801 (2003). 

[4] P. W. Milonni and R. W. Boyd, Phys. Rev. A 69, 023814 
(2004). 

[5] P. W. Milonni, R. Loudon, P. R. Berman, and S. M. 
Barnett, Phys. Rev. A 77 , 043835 (2008). 

[6] D.-W. Wang, A.-J. Li, L.-G. Wang, S.-Y. Zhu, and M. 
Suhail Zubairy, Phys. Rev. A 80 , 063826 (2009). 

[7] F. Intravaia, R. Behunin, P. W. Milonni, G. W. Ford, 
and R. F. O’Gonnell, Phys. Rev. A 84 , 035801 (2011). 

[8] V. Mkrtchian, V. A. Parsegian, R. Podgornik, and W. M. 
Saslow, Phys. Rev. Lett. 93 , 059002 (2004). 

[9] G. Lach, M. DeKieviet, and U. D. Jentschura, Phys. Rev. 
Lett. 108 , 043005 (2012). 

[10] M. S. Tomassone and A. Widom, Phys. Rev. B 56, 4938 
(1997). 

[11] R. Barbieri and J. Sucher, Nucl. Phys. B 134, 155 (1978). 

[12] U. D. Jentschura, J. Phys. A 40, F223 (2007). 

[13] U. D. Jentschura, Phys. Rev. A 81, 012112 (2010). 

[14] U. D. Jentschura and G. H. Keitel, Ann. Phys. (N.Y.) 
310 , 1 (2004). 

[15] Z. C. Yan, J. F. Babb, A. Dalgarno, and G. W. F. Drake, 


Phys. Rev. A 54 , 2824 (1996). 

[16] C. Itzykson and J. B. Zuber, Quantum Field Theory 
(McGraw-Hill, New York, 1980). 

[17] D. P. Craig and T. Thirunamachandran, Molecular 
Quantum Electrodynamics (Dover Publications, Mineola, 
NY, 1984). 

[18] V. Yakhontov, Phys. Rev. Lett. 91, 093001 (2003). 

[19] H. A. Bethe and E. E. Salpeter, Quantum Mechanics of 
One- and Two-Electron Atoms (Springer, Berlin, 1957). 

[20] H. A. Bethe and R. Jackiw, Intermediate Quantum Me¬ 
chanics (Perseus, New York, 1986). 

[21] M. Haas, U. D. Jentschura, C. H. Keitel, N. Kolachevsky, 
M. Herrmann, P. Fendel, M. Fischer, Th. Udem, R. 
Holzwarth, T. W. Hansch, M. O. Scully, and G. S. Agar- 
wal, Phys. Rev. A 73 , 052501 (2006). 

[22] R. Carminati, J.-J. Greffet, C. Henkel, and J. M. 
Vigoureux, Opt. Commun. 261, 368 (2006). 

[23] S. Albaladejo, R. Gomez-Medina, L. S. Froufe-Perez, H. 
Marinchio, R. Carminati, J. F. Torrado, G. Armelles, A. 
Garcia-Martin, and J. J. Saenz, Opt. Express 18, 3556 
( 2010 ). 

[24] A. Manjavacas and F. Javier Garcia de Abajo, Phys. Rev. 
Lett. 105 , 113601 (2010). 

[25] A. Manjavacas and F. Javier Garcia de Abajo, Phys. Rev. 
B 86, 075466 (2012). 

[26] R. Messina, M. Tschikin, S.-A. Biehs, P. Ben-Abdallah, 












10 


P. W. Milonni, and R. W. Boyd, Phys. Rev. B 88, 104307 
(2013). 

[27] G. Lach, M. DeKieviet, and U. D. Jentschura, Cent. Eur. 
J. Phys. 10, 763 (2012). 

[28] T. H. Boyer, Phys. Rev. 182, 1374 (1969). 

[29] J. R. Zurita-Sanchez, J.-J. Greffet, and L. Novotny, Phys. 
Rev. A 69, 022902 (2004). 

[30] G. Lach, M. DeKieviet, and U. D. Jentschura, Phys. Rev. 
A 81, 052507 (2010). 

[31] L. Novotny and B. Hecht, Principles of nano-optics 
(Cambridge University Press, Cambridge, UK, 2012). 

[32] E. A. Power and S. Zienau, Phil. Trans. R. Soc. Lond. A 
251, 427 (1959). 

[33] R. R. Schlicher, W. Becker, J. Bergou, and M. O. 
Scully, in Quantum Electrodynamics and Quantum Op¬ 
tics, edited by A.-O. Barut (Plenum, New York, 1984), 
pp. 405-441. 

[34] W. E. Lamb, Phys. Rev. 85, 259 (1952). 

[35] D. H. Kobe, Phys. Rev. Lett. 40, 538 (1978). 

[36] E. A. Power and T. Thirunamachandran, Am. J. Phys. 
46, 370 (1978). 

[37] P. M. Koch, Phys. Rev. Lett. 41, 99 (1978). 


[38] U. D. Jentschura, J. Evers, M. Haas, and C. H. Keitel, 
Phys. Rev. Lett. 91, 253601 (2003). 

[39] J. Evers, U. D. Jentschura, and C. H. Keitel, Phys. Rev. 
A 70, 062111 (2004). 

[40] U. D. Jentschura, J. Evers, and C. H. Keitel, Laser Phys. 
15, 37 (2005). 

[41] U. D. Jentschura, Phys. Rev. A 69, 052118 (2004). 

[42] D. Meschede, W. Jhe, and E. A. Hinds, Phys. Rev. A 41, 
1587 (1990). 

[43] S. Haroche, in Quantum Electrodynamics and Quan¬ 
tum Optics (Les Houches, Session LIII), edited by J. 
Dalibard, J.-M. Raimond, and J. Zinn-Justin (North- 
Holland, Amsterdam, 1992), pp. 767-940. 

[44] U. D. Jentschura, Phys. Rev. A 91, 010502(R) (2015). 

[45] H. R. Reiss, Limitations of gauge invariance, e-print 
arXiv:1302.1212 [quant-ph]. 

[46] A. Messiah, Quantum Mechanics II (North-Holland, Am¬ 
sterdam, 1962). 

[47] U. D. Jentschura, G. Lach, M. De Kieviet, and K. 
Pachucki, Phys. Rev. Lett. 114, 043001 (2015). 

[48] H. A. Bethe, Phys. Rev. 72, 339 (1947). 


